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Hadronic light-by-light scattering contribution to the muon g — 2* 
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We briefly review the current status of the hadronic light-by-light scattering correction to the muon g — 2. Then 
we present our semi-analytical evaluation of the pion-pole contribution, using a description of the pion-photon- 
photon form factor based on large- Nc and short-distance properties of QCD. We also sketch an effective field 
theory approach to hadronic light-by-light scattering. In view of several still unsolved problems, our conservative 



estimate for the full hadronic light- by- light scattering contribution is a p 



+8 (4) x 10" 



1. Introduction 

The present picture of hadronic light-by-light 
scattering is shown in Fig. 1 and the corre- 
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Figure 1. The hadronic light-by- light scattering 
contribution to the muon g — 2. 

sponding contributions to are listed in Table 1, 
taking into account the corrections made in the 
two full evaluations [1,2], after we had discovered 
the sign error in the pion-pole contribution [3,4]. 
There are three classes of contributions to the 
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Table 1 

Contributions to a M (xl0 10 ) according to Fig. 1. 
The last column gives the result when no form 
factors are used in the couplings to the photons. 



Type 


Ref. [1] 


Ref. [2] 


Ref. [3] 




(b) 


-0.5(0.8) 


-1.9(1.3) 




-4.5 


(c) 


8.3(0.6) 


8.5(1.3) 


8.3(1.2) 


+oo 


/0)Ol 


0.174 a 


-0.4(0.3) 






(d) 


1.0(1.1) 


2.1(0.3) 




~ 6 


Total 


9.0(1.5) 


8.3(3.2) 


8(4) b 





a Only a i exchange. 

b Our estimate, using Refs. [1-3]. 



hadronic four-point function [Fig. 1(a)], which 
can be understood from an effective held the- 
ory (EFT) analysis of hadronic light-by- light scat- 
tering [5,4]: (1) a charged pion loop [Fig. 1(b)], 
where the coupling to photons is dressed by some 
form factor (p-meson exchange, e.g. via vector 
meson dominance (VMD)), (2) the pseudoscalar 
pole diagrams [Fig. 1(c)] together with the ex- 
change of heavier resonances (/o, oi, • • ■) and, fi- 
nally, (3) the irreducible part of the four-point 
function which was modeled in Refs. [1,2] by a 
constituent quark loop dressed again with VMD 
form factors [Fig. 1(d)]. The latter can be viewed 
as a local contribution tjja^ipF^ to g — 2. The 
two groups [1,2] used similar, but not identical 
models which explains the slightly different re- 
sults for the dressed charged pion and the dressed 
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constituent quark loop, although their sum seems 
to cancel to a large extent and the final result is 
essentially given by the pseudoscalar exchange di- 
agrams. Since the models used in Refs. [1,2] are 
not fully consistent with QCD, we take the dif- 
ference of the results as indication of the error 
coming from the model dependence. 

On the other hand, we will show in Section 2 
that the pseudoscalar contribution now seems un- 
der control, due to our semi-analytical calcula- 
tion [3] , based on a pion-photon-photon form fac- 
tor J 7 7r o 7 * 7 » which fulfills the relevant QCD short- 
distance constraints [6], in contrast to the form 
factors used in Refs. [1,2]. These findings are also 
corroborated by an EFT and large- Nc analysis of 
fl LbyL;h a d which allows to calculate the leading 
and next-to- leading logarithms (Sec. 3). 

2. Pion-pole contribution 

The contribution from the neutral pion inter- 
mediate state is given by a two-loop integral 
that involves the convolution of two pion-photon- 
photon transition form factors, see Fig. 1(c). 
Since no data on the doubly off-shell form fac- 
tor ^ r 7r o 7 * 7 » [q 2 , q 2 ) is available, one has to resort 
to models. In order to proceed with the ana- 
lytical evaluation of the two-loop integrals, we 
considered a certain class of form factors which 
includes the ones based on large-iVc QCD that 
we studied in Ref. [6]. For comparison, we have 
also used a vector meson dominance (VMD) and 
a constant form factor, derived from the Wess- 
Zumino-Wittcn (WZW) term. For all these form 
factors we could perform all angular integrations 
in the two-loop integrals analytically [3] . 

In large-A^c QCD, the pion-photon-photon 
form factor is described by a sum over an infinite 
set of narrow vector resonances, involving arbi- 
trary couplings, although there are constraints 
at long and short distances. The normaliza- 
tion is given by the WZW term, J r 7r o 7 * 7 , (0, 0) = 
-N c /(12w 2 F 7r ), whereas the OPE tells us that 



lim J r Tr o 7 * 7 ,(A 2 g 2 , (p — Xq) 



A — >oo 
3 q 2 



T1 
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A2 



1 q-p 



(1) 



A 3 q 2 

In the following, we consider the form factors that 



are obtained by truncation of the infinite sum in 
large- Nc QCD to one (lowest meson dominance, 
LMD), and two (LMD+V), vector resonances per 
channel, respectively: 



■r-LMD / 2 2\ _ - 1 7T 

J- T T o 1 * 1 *(q 1 ,q 2 ) — — 



q\ 



c v 



3 {q\ 



M 2 ){q 2 -M 2 Y 



(2) 



+h 2 qjq 2 + h 5 (qj+qi) + h^j j ((qj - M 2 J 

X (qf - M 2 2 ) {q\ M Vl ) (q 2 - M 2 2 )) j , (3) 

with the constants c v = N c M^/(4tt 2 F 2 ) and 
h 7 = ~N c M^M^/(4:n 2 F 2 ). The parameters 
hi,li2, and h§ in the LMD+V form factor are 
not fixed by the normalization and the leading 
term in the OPE. We have determined these co- 
efficients phcnomcnologically [6,3]. In particular, 
J r 7r o 7 * 7 * (— Q 2 ,0) with one photon on-shell behaves 
like 1/Q 2 for large spacelike momenta, Q 2 = —q 2 - 
Whereas the LMD form factor does not have 
such a behavior, it can be reproduced with the 
LMD+V ansatz, provided that hi — 0. A fit to the 
data yields moreover h$ = 6.93 ± 0.26 GeV 4 . An- 
alyzing the experimental data for the decay tt° — ► 
e+e - leads to the loose bound |/i2|^20 GeV 2 . 

Note that the usual VMD form factor 
^ r (lll 2 2) - l/[(? 2 - M 2 )(q 2 - My)] does 
not correctly reproduce the OPE in Eq. (1). 

After performing the angular integrations, the 
pion-exchange contribution to a M can be writ- 
ten as a two-dimensional integral representation, 
where the integration runs over the moduli of the 
Euclidean momenta 



LbyL;-n 



/>oo />oo 

/ dQi / dQ 2 
Jo Jo 

x£>i(Qi,Q 2 ) fi(Qi,Q 2 ), (4) 



with universal [for the above class of form factors] 
weight functions Wi (rational functions, square 
roots and logarithms) [3] . The dependence on the 
form factors resides in /j. In this way we could 
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separate the generic features of the pion-pole con- 
tribution from the model dependence. This is not 
possible anymore in the final analytical result (as 
a series expansion) for a^ byLi7r ° in Ref. [7]. One has 
to keep in mind that there is an intrinsic uncer- 
tainty in the form factor of 10 — 30 %, furthermore 
the VMD form factor used in that reference has 
the wrong high-energy behavior. 

The weight functions Wi in the main contribu- 
tion are positive and peaked around momenta of 
the order of 0.5 GeV. There is, however, a tail 
in one of these functions, which produces for the 
constant WZW form factor a divergence of the 
form ln 2 A for some UV-cutoff A. Other weight 
functions have positive and negative contribu- 
tions in the low-energy region, which lead to a 
strong cancellation in the corresponding integrals. 

In Table 2 we present the numerical results for 
the different form factors. All form factors lead to 
very similar results (apart from WZW). Judging 
from the shape of the weight functions described 
above, it seems more important to correctly re- 
produce the slope of the form factor at the origin 
and the available data at intermediate energies. 
On the other hand, the asymptotic behavior at 
large Qi seems not very relevant. The results 
for the LMD+V form factor are rather stable un- 
der the variation of the parameters, except for 
h 2 - If all other parameters are kept fixed, our 
result changes in the range \h 2 \ < 20 GeV 2 by 
±0.9 x 10~ 10 from the value for h 2 = 0. 

Thus, with the LMD+V form factor, we get 

a bbyL; "° = +5.8 (1.0) x 10~ 10 , (5) 



Table 2 

Results for a Lb y L ^° for the different form factors. 
In the WZW model we used a cutoff of 1 GeV in 
the divergent contribution. 



Form factor 


a L byL +> x 1Q 10 


WZW 


12.2 


VMD 


5.6 


LMD 


7.3 


LMD+V {h 2 = GeV 2 ) 


5.8 



where the error includes the variation of the pa- 
rameters and the intrinsic model dependence. A 
similar short-distance analysis in the framework 
of large- Nc QCD and including quark mass cor- 
rections for the form factors for the rj and rj' was 
beyond the scope of Ref. [3]. We therefore used 
VMD form factors fitted to the available data for 
J r 7r o 7 , 7 »(— <5 2 ,0) to obtain our final estimate 

LbyL;PS LbyL;jr° , LbyL;r, I . LbyL;T)' I 

U fl — U fl ' U ^ |VMD T U,^ |VMD 

= +8.3 (1.2) x 10" 10 . (6) 

An error of 15 % for the pseudoscalar pole contri- 
bution seems reasonable, since we impose many 
theoretical constraints from long and short dis- 
tances on the form factors. Furthermore, we use 
experimental information whenever available. 

3. EFT approach to a J; byL;had 

In Ref. [4] we discussed an EFT approach to 
hadronic light- by- light scattering based on an 
effective Lagrangian that describes the physics 
of the Standard Model well below 1 GeV. It 
includes photons, light leptons, and the pseu- 
doscalar mesons and obeys chiral symmetry and 
U(l) gauge invariancc. 

The leading contribution to aJ; byLihad , of order 
0(p 6 ), is given by a (finite) loop of charged pions 
with point-like electromagnetic vertices. Since 
this contribution involves a loop of hadrons, it 
is subleading in the large- Nc expansion. 

At order p 8 and at leading order in Nc, we 
encounter the divergent pion-pole contribution, 
diagrams (a) and (b) of Fig. 2, involving two 
WZW vertices. The diagram (c) is actually fi- 
nite. The divergences of the triangular subgraphs 
in the diagrams (a) and (b) are removed by in- 
serting the counterterm \ from the Lagrangian 

£(6) = ( a 2 /47r 2 F ) \ ^75^ ^tt° H , see the 

one- loop diagrams (d) and (c). Finally, there is 
an overall divergence of the two-loop diagrams (a) 
and (b) that is removed by a local counterterm, 
diagram (f). Since the EFT involves such a local 
contribution, we will not be able to give a precise 
numerical prediction for a LbyL;had 

Nevertheless, it is interesting to consider the 
leading and next-to-leading logarithms that are in 
addition enhanced by a factor Nc and which can 
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This cancellation is now also visible in the re- 
vised version of Ref. [8]. In that paper the re- 
maining parts of ci have been calculated: c\ = 
-2x(/i )/3 + 0.237 = — 0.93±^||, with our con- 



ventions for x an d x(M p ) C x P = 1-75_ 1 



4. Conclusions 
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Figure 2. The graphs contributing to a L 1 byL; ' r at 
lowest order in the effective field theory. 



be calculated using the renormalization group [4] . 
The EFT and large- Nc analysis tells us that 



LbyL;had 



(-) J f I m?r± 
\ttJ ] V m,,. 



m K ± 



m. 




m% N, 



c 



16tt 2 F2 3 

2 \ 



, 2 Mo . , Mo . 
m h ci In h Co 

rriu, m u 



tO (-£ x log's) + 0[^N l.,?.-s) 



(7) 



where f(m 7r ±/m fJj ,mx±/m ll ) =—0.038 represents 
the charged pion and kaon-loop that is formally 
of order one in the chiral and Nc counting and /xo 
denotes some hadronic scale, e.g. M p . The coeffi- 
cient C of the log-square term in the second line is 
universal and of order Nc, since F 7r = 0(y / Nc)- 

Unfortunately, although the logarithm is size- 
able, in a^ byL; " there occurs a cancellation be- 
tween the log-square and the log-term. If we fit 
our result for the VMD form factor for large M p 
to an expression as given in Eq. (7) , we obtain 



= (-Yc 



7T/ 

n 



ln2 Mp 



ci In — - 

m. 



Co 



C [3.94- 3.30+ 1.08] 



= [12.3-10.3 + 3.4] x 10 
= 5.4 x 10~ 10 , 
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(8) 



which is confirmed by the analytical result of 
Ref. [7] (setting for simplicity m^o — m^): 



S = [12 - 8.0 + 1.7] x 10~ 10 = 5.7 x 10 



-10 



The pseudoscalar pole contribution aJ; byL;PS 
seems to be under control at the 15 % level. More- 
over, the EFT and large- Nc analysis shows a sys- 
tematic approach to a LtyL;had an( j yi e ^ s i eac ]_ 

ing and next-to-leading logarithmic terms, en- 
hanced by a factor Nc , although these terms tend 
to cancel each other. However, there remains the 
issue of the other contributions, i.e. the dressed 
charged pion and the dressed constituent quark 
loop, see Fig. 1 and Table 1, where model cal- 
culations lead to slightly different results. Taking 
these uncertainties into account by adding the er- 
rors linearly, my estimate for the full hadronic 
light- by- light scattering contribution reads: 

a LbyL;had = + g ^ y JQ-10 _ (9) 
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